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Abstract. Higher order normalizations are performed in the generalized pho- 
togravitational restricted three body problem with Poynting-Robertson drag. 
In this problem we have taken bigger primary as a source of radiation and 
smaller primary as an oblate spheroid. Whittaker method is used to transform 
the second order part of the Hamiltonian into the normal form. We have also 
performed Birkhoff's normalization of the Hamiltonian. For this we have uti- 
lized Henrard's method and expanded the coordinates of the infinitesimal body 
in double D'Alembert series. We have found the values of first and second order 
components. They are affected by radiation pressure, oblatcncss and P-R drag. 
Finally we obtained the third order part of the Hamiltonian zero. 
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1. Introduction 

The restricted three body problem describes the motion of an infinitesimal mass moving 
under the gravitational effect of the two finite masses, called primaries, which move in 
circular orbits around their center of mass on account of their mutual attraction and the 
infinitesimal mass not influencing the motion of the primaries. The classical restricted 
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three body problem is generalized to include the force of radiation pressure, the Poynting- 
Robertson effect and oblateness effect. 



J. H. Poynting(1903) considered the effect of the absorption and subsequent re- 
emission of sunlight by small isolated particles in the solar system. His work was later 
modified by H. P. Robertson (f 937) who used precise relativistic treatments of the first 
order in the ratio of the velocity of the particle to that of light. Chernikov Yu. A. (1970) 
and Schuerman(f980) who discussed the position as well as the stability of the Lagrangian 
equilibrium points when radiation pressure, P-R drag force are included. Murray C. D. 
(1994) systematically discussed the dynamical effect of general drag in the planar circu- 
lar restricted three body problem, Liou J.C.et a/. (1995) examined the effect of radiation 
pressure, P-R drag and solar wind drag in the restricted three body problem. 



Moser's conditions(1962), Arnold's theorem(1961) and Liapunov's theorem (1956) 
played a significant role in deciding the nonlinear stability of an equilibrium point. Moser 
gave some modifications in Arnold's theorem. Then Deprit and Deprit(1967) investi- 
gated the nonlinear stability of triangular points by applying Moser's modified version of 
Arnold's theorem(1961). Maciejewski and Gozdziewski(1991) described the normaliza- 
tion algorithms of Hamiltonian near an equilibrium point. Niedziclska(1994) investigated 
the nonlinear stability of the libration points in the photogravitational restricted three 
body problem. Mishra P. and fshwar B.(1995) studied second order normalization in the 
generalized restricted problem of three bodies, smaller primary being an oblate spheroid, 
fshwar B.(1997) studied nonlinear stability in the generalized restricted three body prob- 
lem. 



In this paper higher order normalizations are performed in the generalized photogravi- 
tational restricted three body problem with Poynting-Robertson drag. Whittaker method 
is used to transform the second order part of the Hamiltonian into the normal form. We 
have performed Birkhoff's normalization of the Hamiltonian. For this we have utilized 
Henrard's method and expanded the coordinates of the third body in double D'Alembert 
series. We have found the values of first and second order components. The second order 
components are obtained as solutions of the two partial differential equations. We have 
employed the first condition of KAM theorem in solving these equations. The first and 
second order components are affected by radiation pressure, oblateness and P-R drag. 
Finally we obtained the third order part H 3 of the Hamiltonian in I^ 2 , I^ 2 zero. 
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Equations of motion are 



dx r\ 
dUi WiN 2 



x-2ny = U Xl where, U x = — — (1) 



y + 2nx = U yi Uy = —L-^-± ( 2 ) 

n 2 (x 2 +y 2 ) (1 -fi) qi n pA 2 

u i = o 1 1 1" ttt y 6 ) 

2 n r 2 2r 2 i 



(x + M^+y 2 , r| = + l) 2 + y 2 , n 2 = l + ^4 



- -9 i 

2 2 



JVi= (» + M)[(» + M)» + w] +i ._ By) ^ a = jd(£±4£±M + ^ + B(3 . + ^ 



= ~ ^g^ 1 M = mi " 2 „ l2 < |j mi,m 2 be the masses of the primaries, A 2 = r " 5r l p 
be the oblateness coefficient, r e and r p be the equatorial and polar radii respectively, r 

/ F \ 

be the distance between primaries, q = (1 — -pr 1 ) be the mass reduction factor expressed 
in terms of the particle's radius a, density p and radiation pressure efficiency factor \ 
(in the C.G.S. system) i.e., q = 1 — 5,6 — -. Assumption g = constant is equivalent to 
neglecting fluctuation in the beam of solar radiation and the effect of solar radiation, the 
effect of the planet's shadow, obviously q < 1. Triangular equilibrium points are given 
by U x = 0, Uy = 0, z = 0, y ^ 0, then we have 



I 1 3/z(l - /i)j/ a;o 2 .T | ^ 



, n^ 2 [2 M -l-Ml^^)f + 7(1-^)4*] ^ 2 (l-f)A 2 | 1/2 

y* = yoii L 5-7; — Vt^ — 2^ — r ( 5 

3/z(l - /x)yg y z 1 



where xo = \ — u, y = ±5(l — ^-) 1 ^ 2 and S = , as in Kushvah & Ishwar(2006) 
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3. Normalization of H 2 

We used Whittaker (1965) method for the transformation of H 2 into normal form. The 
Lagrangian function of the problem can be written as 



1/.2 , -2x , / ■ • n , 71 i 2 , 2x , (i-/")^ . M , M: 



L = -(^ + + n(xy - ±y) + +y a ) + + £ + ^ (6) 

+ ^{ (X + ^ + ^ -narctan. » * 



2r? (x + fj,) . 



and the Hamiltonian is H = —L + p x x + p y y, where p x ,p y are the momenta coordinates 
given by 

dL . W lf dL . VKi 



For simplicity we suppose qi = 1 — e, with |e| << 1 then coordinates of triangular 
equilibrium points can be written in the form 

^ 6 ^ 2 ■ ^ (9 + 7) n Wl -i2L nWl (7) 



2 3 2 3 6V3 27V3 

V3r 2e A 2 2A 2 e (1 + 7 ) 4 7 e -i 

y = < 1 h ^-nWi =nWi > 

2 I 9 3 9 9a/3 27a/3 J 



(8) 



where 7=1 — 2/x. We shift the origin to L4. For that, we change x — > a;* + x and 
2/ — > y* + y. Let a = x* + /x, 6 = y* so that 

Expanding L in power series of x and y, we get 

L = L + L 1 + L 2 + L 3 + --- (11) 
H = H + H 1 +H 2 + H 3 + --- = -L + p x x+p y y (12) 

where L , L\,L 2 , L 3 . . . are 

3 2e 7 e 3 7 A 2 3A 2 e 

Ln = jAo 

u 2 3 3 4 2 ' 

V3nVFi 2 7 nWi neVFi 23e 7 nV^i Tlr 6 

1 = = = nVviarctan- 13 

4 3V3 3V3 54V3 a 
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. f V3 e 5A 2 7eA 2 4nVKi 7 nWi 
~ T" + 3^5 ~ 8Vl + + "9 18~ 

. f 1 e A 2 7nWi 2enVKi 1 

+ y \ 2 ~ 3 ~ T + T2" ~ "HTT + ITT j 

_1 7 9^2 I57A2 35eA 2 29 7 eyl2 3V3nVKi 5enVFi 7 7 ereWi 
~ J <: ~2 + 2 + ~ + ~~ 8 12 12 - + 8 12V3 ~ 4^3 



f 15^42 + V|v42 _ _ 2V - 3i£A2 _ ^ 



43enW 1 
36 



23 7 enWi 1 
36 J 



(14) 



(X ^ ) + n(a# - iy) + y(^ 2 + V 2 ) - Ex 2 - Fy 2 - Gxy 
L 3 = ~ {x 3 T! + 3x 2 yT 2 + 3xy 2 T 3 + y 3 T 4 + 6T 5 } 



where 



E 



= ^{2-6-3^2- 



31A 2£ (69 + 7)^^ 2(307+ 757)^ 



6V3 



27V3 
2(226 + 99 7 )e 

27^ 



(15) 
(16) 



nWi}} (17) 



16 I 



1* 

293A 2 e (187 + 27 7 ) 



6V3 
4(247 + 3 7 )e 



27V3 



+ 7 {6e - — I- v ' ' ' ' nWi = // 

18 6V3 27V3 



W 1 }} 



(18) 



G 



2V3 



nW 1 



2(79 - 77)e 
27V3 



33A 2 e (H-7) TT , (186 -7)e Tlr , 1 
7 {6 - - + 13A 2 - — ^- + '—^nW, - - jj* nW^j 



2V3 



9V3 



(19) 



Ti = 



3 

16 



16 „, 979 , (143 + 9 7 ) „, (459 + 3767) „, 



L„ 4e „ „ 1507 „ (215 + 297) 
+7 14 + T + 25^ - - 



2(1174+ 1697) 
27\ 3 



(20) 
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T 2 



3\/3 
16 



,,16 A 2 367 , 115 1 + 7 „, 959 - 136 7 .„ 

14 - — e + — - - — A 2 e H ' nWi - =— ^nWie 

3 3 18 18V3 27V3 



32e tn , 382 „ (511 + 53 7 ) TTr 

+7 _ + 40A 2 -— ^ + L__J2 nW , 1 



(2519 - 24 7 ) 
27\/3 



nWie 



(21) 



T 3 - 



-9 
T6 



8 203^1 2 625 . (105 + 15 7 ) _ _ _ (403 - 114 7 ) 

r + — - ia a * - V nfl " Siva " Wic 



„ 4e 55A 2 



797 , 197 + 23 7 „, 

— A 2 e + - ^-^nW x 

54 18V3 



(211 - 32 7 ) 
81 V3 



nW\e 



(22) 



-9V3 
16 



8 23.4, 



2 -r 



UA 2 e 



(37 + 7) nWi _(219±M37) nWie 



18^3 



81V3 



88 , (241 + 45 7 ) _ _ _ 



(1558- 126 7 ) 
81 V3 



nW\e 



(23) 



T B = 



2{a 2 + b 2 f 
- 2 (xx + yy ) (ax + by ) (a 2 + b 2 ) 



(ax + by) {3(ax + by) - (bx - ay) 2 } 



(24) 



The second order part H 2 of the corresponding Hamiltonian takes the form 

H 2 = ^±^L + n (yp x - xp y ) + Ex 2 + Fy 2 + Gxy (25) 



To investigate the stability of the motion, as in Whittaker(1965), we consider the following 
set of linear equations in the variables x,y: 
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X 



x 

y 

Px 



and A 



2E G A -n 

G 2F n A 

-A n 1 

—n —A 1 



(27) 



Clearly \A\ = 0, implies that the characteristic equation corresponding to Hamiltonian 
H 2 is given by 



A + 2{E + F + n z )Y + AEF — G+n— 2n 2 (E + F) = 



(28) 



This is characteristic equation whose discriminant is 



D = 4(£ + F + n 2 ) 2 - ±{AEF - G 2 + n 4 - 2n 2 (E + F)} 



(29) 



Stability is assured only when D > 0. i.e 

fj, < ^0 -0.2218959162773076696+2.1038871010983331^2 
+ 0.493433373141671349e^ 2 + 0.704139054372097028nV^i 
+ 0.401154273957540929neT4^i 



where ^ Co = 0.0385208965045513718. When D > the roots ±iu 1 and ±iui 2 (w 1 ,w 2 
being the long/short -periodic frequencies) are related to each other as 



2 , 2 



2 2 
^1^2 



, 76 37A2 83eA 2 2997eyl2 nW 1 5jnWi 53enVKi 
~Y + ~2~ + 12 + 144 24V3 + 8V3 54^3 
h^ 2 nWi \73~ienWi Z^enWx 
24V3 54\/3 36\/3 

27 27 7 2 9e 9 7 e 3 7 2 e 117 7 A 2 2416^ 25157eA 2 



16 16 8 8 8 16 32 

35nW 1 55V3 7 nT4^i 5V3"/ 2 nWi 1277enV^i 



+ ■ 



192 



+ 



16V3 



16 



288V3 



+ 



502l7enT4^i 991j 2 enW : 



+ 



288^ 48V3 

(0 < lu 2 < -^= < wi < 1) 
v2 



(31) 
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From ( |3T)| and (EH) it may be noted that u>j(j = 1,2) satisfy 



1 



4e 107eA 2 2 7 e 1579 7 e,4 2 25nT^i 55 7 nl^i 3809enPFi 



9 27 3 324 

496l7enWi / 16 32e 8 7 e 

486^3 + V 27 + 243 + ~27 ~~ 81 27 
296nW! 10 7 nWi 15892enWi 1864 7 enWi 



27\/3 9V3 486V3 
208^2 8 7 A 2 4868eA 2 563 7 eA 2 



729 



243 



243^ 27 V3 
16 32e 208A 2 
27 



2187V3 729\/3 
1880eA 2 272071^1 



243 



81 



729 



2187V3 



49552enVFi 80 7 enVFi\ 



6561V3 2187%/3 / 

Alternatively, it can also be seen that if u = u)\U)2, then equation 

2 4e 107e^ 2 25nVFi 3809enW^i 

- = H H - 

9 27 27V3 486V3 

55 7 nVFi 4961 7 enWV 



(32) 



1J) gives 



7 



+ 7 



2e 

y " 

16 
27 



27 
1579eA 2 

324 ~9V3 ' 486V3 
32e 208A 2 1880e.4 2 , 320/1^1 15856enWi 
" 243 ^ 



81 



729 



243V3 2187V3 



(33) 



Following the method for reducing iJ 2 into the normal form, as in Whittaker(1965), use 
the transformation 



X = JT where X = 



P % = {2I l oj l ) 1/2 cos 



X 

y 

Px 

Py . 

Q, 



, J — [Jij]l<i,j<4,i T — 



Pi 

Pi 



( — ) 1/2 f 



(i = l,2) 



(34) 



(35) 



The transformation changes the second order part of the Hamiltonian into the normal 
form 

H 2 = cJih - u 2 h (36) 
The general solution of the corresponding equations of motion are 

Ii = const., <j>i — ±uji + const, (i — 1, 2) (37) 

If the oscillations about L4 are exactly linear, the Eq.( [37f represent the integrals of 
motion and the corresponding orbits is given by 

X = J 13 y / 2uj 1 I 1 COS 01 + J14 ^/2uj 2 h COS 02 (38) 



2h 

V = -hi\ sm0i + J 2 

V U>1 

+ J24 \f2L~2UJ2 COS (j) 2 



■ sin 4>2 + J23 v 2/i wi cos < 



(39) 
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where 



Jl3 = 



+^ 

1 

~2fc? 

7 

e 



1 - 



1 

27? 



45A 2 717A 2 e (67+ 19 7 ) „, (431 - 3 7 ) T „ 
e H ^ ;T7 +- + - - ^ ; nVKi - -^nWie 



36 



12V3 



27V3 



3e 

e 
2 



29A 2 (187+27 7 ) nWi _ 2(247 + 3 3 nWie 



36 



3A 2 



12V3 



27V3 



3A 2 



73A 2 e (1-97) „, 53 - 39 7 T „ 

- -— - - + j^nWi + - Ruffle 

24 24V3 

299A 2 e (6 - 5 7 ) 



72 



12V3 



\l\k\ 



3A 2 (33 + 147) w 

— I t= — nVKi 

4 12V3 



54 V3 

w (266 - 93 7 ) 
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54V3 

347A 2 (43 - 87) 



7e 



36 



4V3 



nWi 



(40) 



Jl4 — 



h 


{ 


2uj 2 k 2 




7 

211 


3( 


1 


c 


2fc 2 


2 


+ ^ 
2k 2 




c 


e 




4/|fc 


2 
2 



45A 2 717^ (67+197) „. (431 - 3 7 ) w 
e H — 1 y= — nWi = — nWie 



293A 2 
36 



+ 



2 36 
(187 + 27 7 ) 



12V3 



nWi 



12\/3 
2(247 + 3 7 ) 
27V3 



27V3 



nW x t 



73A 2 e (I-97) „, (53 - 39 7 ) „, 
- 3A 2 — =- + p+nl+i + ^^nWie 



3A 2 



24 
299A 2 e 



24^ 

(6 - 5 7 ) 



72 12^3 
33^2 (1643 - 937) 



nWx 



216V3 



-nWi 



+ 



7e 

4?2fc2 



54 V3 
(268 - 9 7 ) 
54 V3 
737A 2 



72 



nWie 

(13 + 2 7 ) 
73 



nWi 



>(41) 



J21 



4nwi 
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45A 2 
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'2/2 

1 

'2kf 

7 
4fc2 

e 

_ 8lffc2 



3e 



1 + 27f e + 
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36 
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12V3 
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24 
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36 
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12V3 
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27V3 



nWi 



(413 - 3 7 ) 
27^3 



nWie 
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24^ 
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97) TT7" (53 
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54 V3 
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24 V3 
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4V3 



nl+i 
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2l\ [ 36 12^3 
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12V3 27V3 
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7 

e 
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_ 12VT 
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-7 
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2^ 
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+ ^ [51A 2 + ( 14+ ^ n^ 
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1_ 
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7 

2fc 2 

7e 
~8Z 2 fc 2 



3V3 



6^ 



6e + 135A 2 - ^£ _ fci^W, _ g+M nWie 



3e - 18A 2 



9 

55^1 2 e 



2V3 

(1 - 9 7 ) 
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4V3 
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12V3 
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2V3 
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J24 



1 2e + 6A 2 + — — 



4w 2 / 2 /c 2 



(13 + 7) „, 2(79 - 77) „, 



2V3 



-7 



6 + 



2c 



33A 2 e (II-7) 
13A 2 + - 7 J^-nW 1 



9^ 
(186 -7) 



2V3 



211 
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2k\ 



51A 2 



(M + 87) 
3V3 



k\ 



3A 2 



9V3 
(W + 67) 



nW\t 



6V3 



nW 1 



,nr. 808A 2 e (67 +19 7 ) (755 + 19 7 ) 
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2V3 
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3e - 18A 2 -r- - p^nlfi + =-^nVKie 
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12V3 



7e 
21.2 



4Z 2 fc 



2 2V3 



(45) 



with / 2 = 4w 2 + 9, (j = 1, 2) and kf = 2u\ -\,k\ = -2w% + 1. 
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4. Second Order Normalization 



In order to perform Birkhoff 's normalization, we use Henrard's method[Deprit and Deprit 
Bartholome (1967)] for which the coordinates (x,y) of infinitesimal body, to be expanded 
in double D'Alembert series x = J2 n >i^n°> V = Sn>i-^n' 1 where the homogeneous 
components B^'° and B^ 1 of degree n are of the form 

^ Il^Ii ^2 C n-m,m,p,q COS (p<f>i + q<f> 2 ) + S n - m , m ,p,q sin Wl + #2) (46) 
0<m<n (p-.q) 

The conditions in double summation are (i) p runs over those integers in the interval 
< p < n — m that have the same parity as n — m (ii) q runs over those integers 
in the interval — m < q < m that have the same parity as m. Here I\, I2 are the 
action momenta coordinates which are to be taken as constants of integer, </>i, <p2 are 
angle coordinates to be determined as linear functions of time in such a way that <f>\ — 
W i + En>i f2n{h,h), 4>2 = — + J2 n > i 92n (h , h ) where lui , uj 2 are the basic frequencies, 
fm and <72n are of the form 

fin = E / ' 2{n-m).2m,I\ '"^2™ ( 4 7) 
0<m<n 

ff2n = ff , 2(n-m),2m-'l "^2™ ( 48 ) 

0<m<n 



The first order components Bl'° and B?' 1 are the values of x and y given by ( [38]) (I39P. In 

ji,r 



order to findout the second order components B^' , -B^' 1 we consider Lagrange's equations 



of motion 



d_dL ) _dL =Q ± ( ^L)_9L =Q m 
dt [ dx } dx ' dt ( dy> dy 1 J 

£ - 2ny + (2E - n 2 )x + Gy = + 



i.e. } (50) 

y + 2nx+(2F-n 2 )y + Gx = ^ + 

Since x and y are double D'Alembert series, x^x k (j > 0, fc > 0, j + k > 0) and the time 
derivatives x, y, x, y are also double D'Alembert series. We can write 



E 

n>l n>l n>l 



where x, y, x, y are homogeneous components of degree n in l\^ 2 \ I^ 2 i.e. 



d 



\ " ol,0 _ \ "* 



n>\ n>l 



g B i,o 

lf^ + /2 + /4 + 



as 1 - 

— -^(-w 2 +.g 2 +g 4 + ---) 



(51) 



12 



B.S.Kushvah, J. P. Sharma and B.Ishwar 



We write three components x\, x 2 , X3 of x 



3'2 



x 3 



dB\- 
dB, 



LU 2 



1,0 



dB{- 
d<j> 2 
dB. 



DB, 



1,0 



LOl- 



uj 2 - 



DB 



1.0 



dB 



1 

1,0 



2 
,1,0 



w 2 



902 9^1 
1,0 



1.0 



.92 



902 



1,0 



where 



nR i,o x , 9^'° 8B 
DB 3 ' + h^7 32^:- 



d d 

(701 (702 



Similarly three components xi,x 2 ,X3 of i are 



xi = D 2 B\ fi , x 2 = D 2 Bl fi , x 3 = D 2 B\ fi + 2wi/ 2 



2nl,0 



d 2 B[ 
90? 



(52) 
(53) 

(54) 
(55) 



2nl,0 



— 2(jJ 2 g 2 



d 2 B[ 

901 



In similar manner we can write the components of y, y. Putting the values of x, y, x, y, x 
and y in terms of double D'Alembert series in equation ( I50|) we get 



where 



D +2E —1— -An B 



2(1+ ~A 2 ] D — G> B 



X 2 



2 1 + -A 2 ) D + G> B. 



9L, 



,1,0 



D 



2F - 1 - -An I B^' 1 = K 



9a; 



and I2 



D l,0 D 0.1 

x=B 1 ' ,y=B 1 



dLo 



dy 



D l,0 D 0.1 



(56) 
(57) 



These are two simultaneous partial differential equations in B^° and B 2 ■ We solve these 



equations to find the values of B^° and B^' 1 , from Eq. ( |56"|) and (\E 



0.1 



A 1 A 2 S 2 1 ^ = $2, A 1 A 2 B2' i = where A x = D 2 + uj'f, A 2 = D 2 + lu' 2 (58) 



,0,1 



$2 = {D 2 + 2F- n 2 )X 2 + {2nD - G)Y 2 
* 2 = {2nD + G)X 2 - {D 2 + 2E- n 2 )Y 2 



The Eq.f [58)) can be solved for B\'° and B^ 1 by putting the formula 



(59) 
(60) 



1 



AiA 2 



cos(p0i + #2) 
or 

sin(p0i + g0 2 ) 



1 



cos(p0i + g0 2 ) 



A 



p - q [ sin(p0i + 902 ) 
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where 

A P , g = [u\ - (wip - u 2 q) 2 ] [ul - (wip - w 2 g) 2 ] 

provided A p , q ^ 0. Since Ai^ = 0, A ,i = the terms cos <p\ , sin <p\ , cos 2 , sin 02 are the 
critical terms, $2 and \&2 are free from such terms. By condition(l) of Moser's theorem 
kiui + k 2 u 2 ^ for all pairs (ki,k 2 ) of integers such that \k\\ + \k 2 \ < 4, therefore 
each of uj\ , w 2 , 0J1 ± 2u 2 , u 2 ± 2wi is different from zero and consequently none of the 
divisors A ,o, A , 2 , A 2i0 , Ai^i, Ai__i is zero. The second order components B^.B^' 1 
are as follows: 



B. 



1,0 



and 



B 



0,1 
2 



where 



+ 



rih + r 2 I 2 + r 3 h cos 20! + r A I 2 cos 2<j) 2 + r 5 I^ 2 I^ 2 cos(0i - 2 ) 



r 3 = 



1/2 r l/2 i 



1/2 r l/2 



i>\ + 02) + ^7/1 sin 20i + r$I 2 sin20 2 
sin(0i - 4> 2 ) + r w l{ /2 ll /2 sin(0i + <t> 2 ) 



(61) 



= -{sih 



s 2 I 2 



S3/1 cos 20i + Sil 2 cos 20 2 + s^T^ 2 !^ 2 cos( 



+s 6 i 1 1 / 2 / 2 1 / 2 cos(0i + <j) 2 ) + s 7 Ii sin 20i + s$I 2 sin20 2 
+ s 9 l\ /2 ll /2 sin(0i - 02 ) + si / 1 1/ ' 2 / 2 1/ ' 2 sin(0i + 2 )j 

= ^| I J1V1F4 + J13 J 23 ^l F l + + J 2 2 3 W1 ) i^' 



r 2 



l w 2 



3w 2 (4w 2 - lu 2 ) 



T 2 

J 21 
U)l 

T 2 

J 21 



J14W2-F4 + J14 J 2 4W2-Fi + ( — + jf 4 W2 ) F 4 ' 

LU 2 



8ivfj 21 (J 13 F{ + 2J 23 F{') + Au 2 



(62) 



(63) 



(64) 



- J23W1 \F 1 



J 23 uji )F[' 



{J13F2 + J 23 F 2 )J 13 UJi 
2cJ! J 2 i(Ji 3 F^ + 2J 23 F^) - uoxJuiJ^Fi + J 23 ^')wi 

(65) 



T4 



3uj 2 (Alu 2 - LU 2 ) 



8loIJ 22 (J u F[ + 2J 24 F[') - 4lu 2 



(JuF 2 + J 2 aF 2 )Julu 2 



( 111 


— J 2A L0 2 


K 


V ^2 




( J 22 


- J 24 lo 2 




V ^2 







2lo 2 J 22 (J u F; + 2J 2i F^) - LU 2 J u (J 14 Fi + J 2i F^)uj 2 

(66) 
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1 



WiW2(2wi + W2)(4wi + 2w 2 ) 



(wi + w 2 ) 3 



{Jl 3 J 22 (^) 1/2 -Jl 4 J2l(^) 1/2 }^' 



W 2 s 



W 2 



_ 2 { J 21 J 24 (^)V2 _ J 22 J 23 (^)V2} F « 

^1 <^2 



-(Wl + W 2 ) 2 



{2{ J13 J14F2 



J21J: 



+ (^13^24 + Jl4 JM) h^)' /2 |{ (J^jT/2 + ^3J24^1^2) 1/2 }^' 
~(Wl + ^) [{ J13 J 22 (^) 1/2 - Jl 4 J 2 l(^) 1/2 }^ - 2{ J 21 J 24 (^) 1 / 2 

{2{ J13J14F4 + (J13J24 + Ji4J 23 )i ;, l}Kw 2 ) 1/2 
+ J23 J 2 4(wio; 2 ) 1/2 }Fi / 



J 2 2J 23 (-) 1/2 }F ; 
W 2 



2{ J2lJ22 



(^l^ 2 ) 1 /2 



(67) 



-1 



(wi - uj 2 y 



wiw 2 (2(jJi — w 2 )(4o;i — 2w 2 ) 

+ 2{j 2 lJ24(-) 1/2 + J22J23(-) 1/2 }F 1 " 



U>2 Wl 



{2{ J 13 J 14j F 2 + (J 13 J 24 + Jl4J23)^2}K^2) 1/2 



- J23J 24 (wiW2) 1/2 }i ?, 2 ' 



- (u>i - U) 2 ) 



+ (uj 1 - 0J2) 

„f ^21 ^22 
1 (W1W2)V2 

- Ji 4 J 2 i(^) 1/2 }F 3 ' + 2{ J 21 J22A 1 / 2 + J22 J 23 (^) 1/2 }F3" 

{2{j 13 Jl4i ? 4 + (^13^24 + ^14^23)^}(W1W2) 1/2 



{Jl3*2(^) 1/2 



W 2 ' 



(68) 



3w 2 (4w 2 _ w ^ 
-2wi 



./,.,(./! i r, -- - ( ^ ./ 2 : .~, . ' 



-4o> 2 J 21 (J 13 F 2 + Jasi^V + J2i(Ji 3 Fl + 2J 23 Fi') 



(69) 
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3a>|(4a;| - u){ 

+4^2^22(^14^2 + 2J 2 aF%)lU2 ~ 2iU 2 



— - ^24^2 ) F% 



Jli{JliF\ + J24F[)uJ 2 — ( — — — J| 4 t^2 )F[' 

\ W2 

^2 Jl4(^14-F3 + J24-F3) 
J22{JuF i + 2 J24-F4 ) > 



wiw 2 (2wi + W2)(wi + 2o> 2 ) 



(wi + w 2 ) 3 



{2Ji3 J14F1 



+ (Jl 3 J 2 4 + Jl4J 23 )i ;, l'}(wiW2) 1/2 + 2{ 



-(Wl + cj 2 ) 2 



{Jl 3 J22(^) 1/2 -Jl 4 J2l(^) 1/2 }^ 
U>2 ^1 



J21J22 

/ w 2 < 



+ J2 3 J2 4 K^ 2 ) 1/2 }Fl" 



_ 2 {J 21 J 24( ^)V2_ J 22 J 23( ^)V2}^ 

Wi Co>2 



■(wi + W 2 ) 



{2{ J 13 J 14 F 3 + (J13J24 + Ju^s)^}^!^) 172 



+ 2 { / 2lJ U/ 2 + J 23J2 4 (^2) 1/2 }^l - [{Jl 3 J 22 (^) 1/2 



(CJH02) 1 / 2 



U>2 



~jM^) i/2 }n 2{ j 21 j 24 (^) i/2 - J22 j 23 (^) i/2 k 

Wl Wl W2 
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no 



1 

WlW 2 (2^i — UJ 2 )(2uJ2 — Wl) 



(ui - uj 2 ) 3 



{2J13 J14F1 



+(Ji 3 J 2 4 + JiMFi}^) 1 / 2 - 2{ / 2lJ " 2 1/2 - J2M^ 2 ) 1/2 }F^ 

(UJIUJ2) ' 



-(u>i - UJ 2 ) 



{Jl 3 J 2 2(^) 1/2 -Jl 4 J 2 l(^) 1/2 }^ 



^2 xl/21 



+2{J 21 J 2 4(^) 1/2 + J 22 J 23 (-) 1/2 }^' 



"2{ 



(wi - w 2 ) 
J21J22 



{2{J 13 J U F 3 + (J 13 J 24 + Jl4J 2 3)^3}K w 2 ) 1/2 



(wiO*) 1 / 2 

{Jis^) 1 ' 2 -^!^) 1 / 2 }^ 



w 2 



+2{J 21 J 24 (^) 1/2 -J 22 J 2 3(^) 1/2 }^' 



0J 2 



(72) 



We can write expressions of with the help of replacing by Gj, by G[ and F/' 
by = 1,2,3,4), where 



F-i = 



—nWie 



(73) 



F 2 = — 



16 „ „ 979 , (143 + 9 7 ) (555 + 3767) 

-e + 6A2 - - A 2 e + + '-^^nW.e 



L„ 46 „r . 1507 , (215 + 297) — — - 

+7 14 + T + 25A 2 - — A 2 , - ^-^^r 



2(1174+ I697) 
27V3 



nWie 



(74) 



^3 



3V3 



16 



AA 16 23yl 2 104, 115 1 + 7 2 439 - 687 

14 - ^-e + — — ^^ 2 e H v ' nVKi =-^-nWie 

3 2 9 18V3 27V3 



[32e AnA 310, (511 + 537) (2519- 249 7 ) TTT ] 

+ 7 | - + 40^ 2 - - A 2£ + - ^ ^ lCj 



(75) 
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-3 
256 



™„ 17801A 2 , (2821 + 189 7 ) Tir (23077 + 9592 7 ) TT , 
364 + 420A 2 — -A 2 e + ± = '^nW x - - = -nW x e 



„„ „„ lOOe 849A 2 59 „ 
+28 7 23+ — + — + T A 2C - 



3V3 
(125 + 38 7 ) 
6V3 



27V3 



(87613- 213 7 ) 
27V3 



(76) 



3V3 
~T6~ 



nWie 
3V3 

,, 16 , 1367 , 115(1 + 7) 



(77) 



^ (863- 136^ 
27v 3 



f 32e 382 „ (511 + 53 7 ) Tlr (2519 - 24 7 ) 

+ 7 | - + 40A 2 - -A 2 e + ^-^^i - 2?V 3 7 W 



(78) 



-9 



8 203A 2 721 (105+ 157) „, (319 - 114 7 ) 

s £ + — - n ^ - Svr^ 1 - Siva W l£ 



[ 4e 173A 2 781 A (197 + 23 7 ) 



(265 - 32 7 ) 
81V3 



nW\€ 



(79) 



-3V3 
16 



392 _ 532e 1918^2 _ 28582^ (203 + 1211 7 ) 

3 3 9 9^3 



949 + 4378 7 ) TTr na J 108e 4037^ 2 611. (8397+ 919 7 ) 
27\/3 I 7 84 21 84V3 



(92266- 1869 7 ) 
27V3 



nW\e 



(80) 



nWie 



(81) 



-9 
32 



203A 2 
6 



4e 55A 2 



9 



625 , 105+ 15 7 307- 114 7 w 

— A 2 e - 7 =-^-nWi - =— ^nT+ie 

54 18V3 8lV3 

797 , (197 + 23 7 ) (211 - 32 7 ) 1 



(82) 
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8 55A 2 134 „ 37 + 7) 93 + 226 7 ) 

2 - -e + — - ^-^2£ - T^nWi - =-^nWie 

3 6 3 18V3 81V3 



169, 241 + 45 7 1558 - 126 7 

+7< 4e + — A 2 e + 1 =- LL nW 1 - ^ -= — -nW^ 

I 27 18V3 81V3 



(83) 



9 
256 



212 295(L4 2 1370^ 2 , (771 + 237 7 ) 2(1907 - 984 7 ) 

e H A 2 e = nWi = nW\e 

3 3 27 9V3 81V3 



11 4e 152yl 2 36965 , (2569 + 277 7 

+287^ — + — — ^ - -T7--A 2 e + - =-^nW x 

17 9 7 504 252V3 



, (22603 + 4396 7 ) 



1134 V3 



(84) 



—nWie 



(85) 



_3_ 

32 



,, 16 . 1367 , 115 1 + 7 863- 136 7 

14 - — e + A 2 - — -A 2 e + v -= ' nW 1 - - ^-^nW^ 

3 18 18V3 27\/3 



32e tnA 382 A (511 + 53 7 ) _ _ _ (2519 - 24 7 ) 
+7 - + 40A 2 - -A 2 e + K —^-nW, - ^^ "^ 



(86) 



3V3 



16 fl . 907A 2 A (143 + 9 7 ) jjt (477 + 403 7 ) _ _ _ 



I 4e 71A 2 1489 . (215 + 29 7 ) T „ 



2(1174+ 1697) 
27V3 



(87) 



3V3 
~256~ 



84 + 52, + 212^ 2 - 267A 2£ + ^t^W - ^l+^n^e 



3V3 



27V3 



+7 32, + 156A 2 + 649A 2£ - n Wl + n Wl e 



3V3 



27V3 



(88) 



G[ = 



-nW\e 

~7T 



(89) 
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9 
16 



8 203A 2 625 t 



m + i5,)_ nWi _(m-u^) nWie 



„ 4e 55^2 797 „ 

-7< 2 A 2 e + 

'] 9 2 54 



18^3 
(197 + 23 7 ) 
18V3 



81V3 



(211-327) 
81\/3 



(90) 



3V3 



14 e 

3 



65A 2 
6 



1439 , 115(1 + 7) „, (941- 1187) _ _ _ 



18 



18V3 



27V3 



32e „„„ 310 , (511 + 537) TTr (251- 24 7 ) 



(91) 



-9 
128 



847 A 2 , 2(28 + 7) 4(2210 - 69 7 ) Tir 
12e - 287A 2 + - A 2 e - — — ^-nW x - = — -nWie 



V3 



27V3 



-7<96+^ + 135A 2 



2320 „ (497- 1237) TTr 
^ 2£+ 3V3 nWl 



4(17697 + 32 7 ) 
27V3 



(92) 



G'l = 



—nW\€ 



(93) 



9V3 
"32" 



2-- e+ — -UA 2 e 



(37 + 7)^ i _(123 + M97)^ ie 



18V3 



3V3 



1 , 88A 2 (421 + 45 7 ) „, (1558- 1267) „. 
+l{ 4e + — 7^r~ + - — k nWi - - ^—^ — -nWxe 



27 



18\/3 



81V3 



(94) 



-9 
~L6 



203A 2 589 



54 



A 2 e 



5(51 + 27) 
18V3 



(349 - 282 7 ) 



8lV3 



nW\€ 



L 4e „ 412 „ (197 + 23 7 ) TTr 
+ 7 2- ¥ -26A 2 --A 2 e+L__J2 nWl 



(211 -327) 
81 V3 



nW\e 



(95) 
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-9^3 
256 




(32 7 ) 



3\/3 



nWi 



2(1529 + 450 7 ) 
27V3 



nWie 



+7<8e 



749^2 808 




(109 - 4O7) 
3V3 



nWi + 




(96) 




We can verify that in Eq.( [97)) , A3.0 vanishes independently as in Deprit and Deprit 
Bartholome(1967). Similarly the other coefficients ^2,1, ^4.1,2, A), 3 are also found to be 

1 /2 1 /2 

zero independently. Hence the third order part H3 of the Hamiltonian in I x , I 2 is 
zero. 



Using Whittaker(1965) method we have found that the second order part H2 of the 
Hamiltonian is transformed into the normal form if 2 = wiii — u^i^- The third order part 

1 /2 1 /2 

of the Hamiltonian in Jj/ , J 2 ' is zero. 
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